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Abstract. We study the fluctuations of a stochastic Maxwell-Lorentz particle 
model driven by an external field to determine the extent to which fluctuation 
relations are related to large deviations. Focusing on the total entropy production 
of this model, we show that, although the probability density of this quantity 
globally satisfies (by definition) a fluctuation relation, its negative tail decays 
exponentially with time, whereas its positive tail decays slower than exponentially 
with time because of long collision-free trajectories. This provides an example of 
physical system for which the fluctuation relation does not derive, as commonly 
thought, from a probability density decaying everywhere exponentially with time 
or, in other words, from a probability density having a uniform large deviation 
form. 



PACS numbers: 02.50.-r, 05.10. Gg, 05.40.-a 



1. Introduction 



The fluctuation relation (FR) is an important result of nonequilibrium statistical 
mechanics, expressing a general asymmetry for fluctuations of systems driven in 
nonequilibrium steady states. Such an asymmetry typically applies for observables 
At integrated over a time t and implies that 

P(A t /t = a) _ cta+0 ( t ) /ix 
P(A t /t = -a) ' W 

where P(-) denotes the probability density function (pdf), c is a constant that does not 
depend on t nor a, and o(t) stands for sublinear, i.e., sub-extensive corrections in t. The 
exponential dominance of positive fluctuations over negative ones expressed by this 
result has been widely studied for chaotic and stochastic systems [l]-[5] , as well as for 
different physical observables of these systems, such as the entropy production, particle 
currents, and work- like quantities [6 -12 



see 



13 - 15 for reviews. Observables that 



satisfy the FR of Eq. |l]) have also been measured experimentally, e.g., in manipulated 
Brownian particle experiments and noisy electrical circuits [16 - 20 
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For all of these systems, the exponential form of the FR is known to arise because 
of two fundamental properties of the pdf P(A t /t): (i) it satisfies a large deviation 

i.e., 

lim — \nP(A t /t = a) = 1(a) (2) 



principle (LDP) |2l}|23 
or equivalent ly, 

P(A t /t = a) = e -*n«)+°W ) (3) 
and (ii) the so-called rate function 1(a) has the symmetry property: 

I(-a) - 1(a) = ca. (4) 

These two conditions, with 1(a) different from and oo, are known to be sufficient for 
A t ft to satisfy an FR (see, e.g., [23] and references therein). 

In this paper, we show with a specific physical model that the FR can also arise, 
unexpectedly, from a pdf that does not have a leading exponential scaling in t, as in 
Eq. (J3j) , and so does not satisfy a "standard" or "classical" LDP 1 24 . By considering a 
stochastic Maxwell-Lorentz gas driven by an external field, we show that the entropy 
production calculated over a time t satisfies (by definition) an FR, even though the 
far positive tail of its pdf scales exponentially with \fi rather than t. This provides a 
physical example for which an FR arises not from a single, uniform LDP as above, but 
from different large deviation scales (here two scales) , which can be fully characterized 
only by explicitly calculating the o(t) correction term in Eq. ^. In our model, the 
crossover to the \ft scaling is physically related to long, collision-free accelerations of 
the gas' particles under the external field, which lead to large positive fluctuations of 
the entropy production. 

It is important to note that the features of the model that we study are different 



from those studied in the context of so-called extended FRs 25 and of boundary-term 
effects in FRs 26 -29 . In those studies, the LDP condition of Eq. ^ is satisfied, but 



the associated rate function does not satisfy the symmetry Q . Our results are also not 
related to anomalous FRs, which arise when there is no LDP because the probability 



density of interest has power-law tails 30-32 . Here we study a non-trivial case of 



non-uniform large deviation scaling, which has the interesting feature of preserving 
the exact exponential scaling form of the FR. The relevance of this result for more 
general collisional models is discussed at the end of the paper. 

2. Model and relevant stochastic variables 

The model that we consider is a stochastic Maxwell-Lorentz gas previously analyzed 

in [33 34 , consisting of a probe particle of mass m whose velocity v changes because 

of collisions with particles from a bath and acceleration due to an external force field. 

The collision process is represented by interactions with a bath of scatterers 

of mass M equilibrated at temperature T. Collisions with the scatterers change 

instantaneously the probe's velocity from v to v' according to the rule 

. , . (-a 

v = 71; + (1 - j)V, 7 = 73^, (5) 

where a g [0,1] is the restitution coefficient (a = 1 for elastic collisions), £ = m/M 
is the mass ratio, while V is the velocity of the scatterer, taken to be distributed 
according to the Gaussian pdf: 

Pscatt(^) = X [^e-^ 2 (6) 
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with q = M/(2T). The fact that the pdf above is independent of the probe's velocity 
means physically that the scatterers do not keep any memory of their collisions with 
the probe and, so, that the bath of scatterers thermalizes rapidly, relative to the typical 
interaction time, with the probe particle. 

During a time r between two consecutive collisions, the probe performs a 
deterministic acceleration under the influence of an external field £ . To simplify the 
model, we assume that r is taken for each flight from the exponential pdf, 

P T (r) = ie-^, (7) 

which is independent of the relative velocity of the particles. This is a simplification 
compared to collisional models with hard-core interactions, where the rate of collisions 



is proportional to the relative velocity \v — V\; see, e.g., 35 36 . 

Under the collision rule and acceleration defined above, the evolution of the 
velocity pdf P(v, t) of the probe particle is described by the following linear Boltzmann 
equation: 

T c d t P{v, t) + r c £d v P{v, t) = -P(v, t) + ^— J du P(u, t)P scatt (^yJ , (8) 
where t c appears as the mean collision time. Various properties of this integro- 



diffcrential equation are discussed in 33 34 . In particular, for the case M > m 



and C, = a, i.e., 7 = 0, which implies the simple collision rule v' = V, the stationary 
Boltzmann equation can be solved analytically to find 

'-<•> = VI f* -"- , "- > ' = 5'^'* (^f ) • < 9 > 

with b = 1/(£t c ). This is the case that we consider throughout this paper. For the 
general case 7 ^ 0, an explicit form of the solution is not available, although it the 
can be written as a series expansion in Fourier space [34] . 

It is worth noting that the collisional model described by Eq. (Jsl) belongs to 



a more general class of models recently discussed in 37 , whose collision integral 
includes a term of the form \v — u\ v in the kernel and a scatterers' pdf generalised to 
-PscatttX) ~ e^l^l^ (here we consider v — and \i — 2). In 



37 



it is shown that the 

parameters space (y, /i) is divided by a "transition line" in two regions: fi < v + 1 
corresponds to models where the interaction with the scatterers is strong enough to 
guarantee thermalization, while for fi > v + 1, which is the case considered here, the 
accelerated particle takes a stationary pdf which is far from that of the scatterers. 
The model that we study is also in the same class of the original Lorentz gas |38| and 



the frequently-studied Sinai billiard 39 , for which v — 1 (hard spheres) and /j, = 00 
(scatterers at rest). 



3. Total entropy production 

The quantity that we study for the purpose of comparing the FR and the LDP is 
the total entropy production As to t associated with the velocity v(t). This quantity is 
defined in the standard way as 



(10) 
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where P({v( s )lo) and P({v(s)Y ) are, respectively, the pdf of a path {v(s)}q spanning 
the time interval [0,t] and of the time-reversed path {i>(s)}o = {— v(t— s)}q [5J. From 
this definition, it is easy to check that (see [5|[l5]) 



P(As tot = a) = e a P(As tot = -a), (11) 

so that 

P(As tot /t = s) = u 

P(As tot /t = -s) ' 1 ' 

for all i > and s £ K. Therefore, As tot /i satisfies as announced an FR with c = 1, 
which holds exactly in this case as there is no o(t)-corrections to Eq. ([!]). 

At variance with this result, we show in the next section that P(AstotA = s ) 
itself is not everywhere exponential with t, which means that the FR in this case does 
not derive, following the introduction, from an LDP having a uniform scale or speed 
t, but from a non-trivial combination of LDPs showing different scalings with time. 



To derive these LDPs, we use in the next section a result of 34 showing that 
Astot can be decomposed as 

W 

As tot (t) = — + B (13) 

where 

W(t) = mS v(s)ds = m£^Xi (14) 
Jti i=i 
is the work done by the external field between the first and the last collisions, and 

is a boundary term involving only velocities close to the start and to the end of the 
time interval [0,t]. In the above formulas, n is the (random) number of collisions 
up to time t, 9 — TC, is the energy scale corresponding to equipartition with the 
scatterers (which is not reached because of the external field and the dissipation in 
collisions), while v(tf) and v(t~) are the velocities of the particle before and after the 
z-th collision, respectively. Moreover, 
pti+i 

v{s)ds (16) 

represents the probe's displacement between the times U and ti+%. For 7 = 0, the 
post-collisional velocities v(tf ) are extracted from the Gaussian pdf whereas the 
pre-collisional velocities are distributed according to the pdf 

4. Center LDP for the total entropy production 

Rate functions of random variables are often obtained by calculating their 
corresponding scaled cumulant generating function (SCGF) [23]. For the total entropy 
production, the SCGF is defined by the limit 

A Astot (fc)= lim ^ln( e fcAs -). (17) 
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Following the Gartner-Ellis Theorem 21-23 , it is possible to obtain the rate function 
I(s) of As tot /t, defined by 

I(s) = lim -i]nP(Aa tot /t = «) ) (18) 

t— >oo t 

provided AA Stot (k) is differentiable and steep in the interior of its domain of existencejf] 
Under these conditions, we then have that I(s) is the Legendre-Fcnchel (LF) transform 
of A Astot (fc), i.e., 

I(s)=max{sk-\ A8tot (k)}. (19) 

fcfcK 

The calculation of \As tot (k) is simplified here thanks to the decomposition of 



Eq. (13) and the fact that the work W(t), which is an integral of a random process 
over the time interval [0,t], is asymptotically independent of the boundary term B 1 
which involves only the limits of this integral. Defining the SCGF of the work by 

X w {k) = lim -\n(e kW / e ) (20) 

t—>oo t 

and that of the boundary term by 

X B (k)= lim jln(e kB ), (21) 
we must then have 

A Astot (fc) = X w (k) + X B {k). (22) 



The calculation of Xw(k) is carried out in Appendix A.l with the result 



X w (k)= Gxik) - , (23) 

km£ 

where G x (k) = (e~s~ x ) is the generating function of the displacements Xi appearing 



in the expression of the work, Eq. ( 14 ) . This generating function cannot be obtained 



in closed form; however, it can easily be evaluated numerically, and from its integral 



representation, shown in Eq. (A. 4), we find that Xw(k) < oo for k € (—1,0] and 



Xw{k) = oo otherwise. The domain of existence of Xw{k) is thus (—1,0]. Moreover, 



from Eq. (A.4) we obtain the FR-like symmetry 

X w (k) = X w (-k- 1) (24) 

for all k G K if we accept the equality oo = oo. 

The SCGF of the boundary term B is calculated in a different way by rewriting 
it as the sum B = b± + b n of two asymptotically independent terms: 

6i = S u2(tr)+lnP(w(0)) 

b n = -^v 2 (t+)-lnP(-v(t)). (25) 

Since these terms do not scale extensively with t, their SCGFs must vanish when it 
exists, so that we only need to determine their domain of existence. This is done 
in | Appendix A.2| The result that we obtain is that the SCGF of &i converges for 
k G (—oo,0], while that of bi converges for k G (—1,0]. Combining these domains, 
i.e., taking their intersection, we then find 

**(*) = { ° (26) 
v ' I oo otherwise. 

% The steepness condition means essentially that the derivative of the SCGF diverges at boundary 
points of its domain; see Examples 3.3 and 4.8 of [231 for a complete explanation. 
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Thus As(fc) has the same domain as Xw(k), which means that we can finally write 
\As t Jk) = \ w (k) (27) 



with Xw(k) given again by Eqs. (23) and (A.4). This shows that only the work plays 



a role in the SCGF of the entropy production - the boundary term is irrelevant. 

In Fig. [lja) we plot AA Sto t(^) by numerically computing the integral in Eq. (A.4) 



for several values of the field 8. Notice that, as expected from the symmetry (24), 
AAstotW is symmetric with respect to k = —0.5, and that its right-derivative at 
k = — 1 and left-derivative at k — increase in magnitude as the field £ is increased. 



From Eq. (A.4), we actually find 



A Astot (fc = 0") = -A' Astot (fc = = (28) 

This value is important for our discussion and will be denoted next by s*. 

Having the SCGF of the entropy production, we can now in principle obtain its 
LDP using the Gartner-Ellis Theorem. However, the problem arises that the SCGF 
is not steep, so it is not clear whether the rate function I(s) is the LF transform of 
the SCGF. What is known at least is that this transform yields the correct I(s) for s 
in the image of the derivative of \As tot (k), which here corresponds to (— s* , s*). Over 
this interval, we can indeed express the rate function as the Legendre transform of the 
SCGF §] so that 

I(s) = kX' Ast Jk)-X Astot (k) (29) 

for 

s = X' Ast Jk)e(- S * lS *). (30) 

The result of this Legendre transform is shown in Fig. [ljb). It is worth noticing that 
I(s*) = 0. Moreover, from the symmetry of Eq. p4| , we readily obtain 

I(-s) - I(s) = s, (31) 

so that the original FR symmetry of Eq. Q is satisfied for s e (— s*,s*). Over this 
interval, we have in fact 

P(As tot /t = s) = e- tI ^+< t \ (32) 

so that both the LDP and the FR for the entropy production are exponential in i, in 
agreement with the scenario outlined in the introduction. 

Outside the centered interval (— s*,s*), one cannot be certain that the rate 
function is the LF transform of the SCGF because this transform is one-to-many for 



non-steep functions 23 . Other arguments are given in the next section to obtain the 
rate function outside this interval, with the result that P(Astot/t = s) should decay 
with \ft rather than t for s > s* , implying that I(s) = for s > s* and I(s) — — s for 
s < —s*. This is consistent with the LF transform of AAs tot (fc) and is supported by 
simulation results presented in Fig. [ljb). 

§ See Example 3.3 and Sees. 4.1 and 4.4 of [23] for more details. 




Figure 1. (a) SCGF for the total entropy production, Eq. < |22[ i. for the parameters 
m = 1, M = 2,r c = 1,T = 1,7 = and various values of the field £. 
(b) Numerical estimates of the rate function I(s) obtained at different times 
measured in units of the mean collision time r c . The numerical simulations 
were obtained, following the method described in [34] , with the parameters 
m = 1, M = 2, T = 1, r c = 1, 7 = 0, E = 1. The t = oo black line represents 
the rate function obtained from the Legendre transform of the SCGF; the red 
line shows the part I(s) = —s for s < — s* , while the dotted line shows the part 
I(s) = for s > s*. 



5. Far tail LDPs 

The typical trajectories of the system that give rise to fluctuations of As to t close to 
its typical value s* are those that involve many collisions (of the order of </r c ) over 
the time t. To study the large fluctuations of this quantity away from s* - at least 
positive ones - we must therefore study long ballistic trajectories of the probe particle 
involving few collisions. Physically, these trajectories, exponentially distributed in 
time, have the effect of producing a "slower-than-exponential" or "fat" LDP in the 
far positive tail of the entropy production pdf jjj] 

To see this, consider the extreme case where no collision occurs in [0,t] so that 



t. Then, as shown in 34 , the entropy production is exactly given by 



A,.^.»« (33) 

where v(Q) is distributed according to the stationary pdf of Eq. ^ and v(t) — v(0)+£t. 
From the expression of the stationary pdf, it is easy to see that its negative tail has 
the form of a Gaussian, which implies 

P(-v(t)) - P stat (-v(t) - £t) ~ e -^m+£t)\ (34) 



|| Note that this fat tail does not induce any anomalous diffusion, i.e., one has ([x(t) — (x(t))] 2 ) ~ t 
at large times [34|[40| . 
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for v(t) ^> 1 and £ > 0, so that 

A Stot ~ q£V = ^£ 2 i 2 (35) 

if we retain in the entropy production only the dominant positive term that scales 
with t. The probability of this "no-collision" event is given by the exponential pdf of 
t shown in Eq. Q, so that by change of variables we find 



with subexponential corrections in \fi. 

The same result is obtained if we consider a long ballistic trajectory with one 
collision at either end of [0,t] or two collisions close to the start and end of this time 
interval. In the latter case, the dominant part of the entropy production is the work 
done over the displacement 

x = vt + t 2 £/2 ~ i 2 £/2, (37) 

which leads us to the same asymptotic for Astot a s that shown in (|35[) and, 



consequently, the same LDP shown in (36). Hence the extreme positive fluctuations 



of the entropy production are governed by a "fat" LDP, which scales with y/t rather 
than t. This explains why the SCGF of the entropy production diverges for k > 0: the 



stretched exponential tail of the pdf ( 36 ) is subleading compared to the exponential 
term in the SCGF ( 17), and so cannot compensate for the divergence of this term for 
k > 0. 

These simple arguments cannot be used to determine the precise value s at which 
the change of scaling from t to yi occurs. However, if we assume that there is no 
additional scaling and that the pdf of the entropy production is unimodal, then the 
divergence of AAs tot (k) for k > together with the value s* of its left-derivative at 
k = imply that this crossover value must be s* itself. In this case, if we take the 
large deviation limit with the scaling t, as in the previous sections, we find J(s) = 
for s > s* , as shown in Fig. [ljb). This result is also consistent with the divergence of 
AAstot(^) f° r ^ > 0, though we must stress again that a zero rate function is in general 



only an artifact of taking the large deviation limit with the "wrong" time scaling 23 



it is a signal that a "fat" LDP governs the large fluctuations of the entropy production, 



which for our model has the form shown in ( 36 ) . 

To find the LDP scaling associated with the negative fluctuations of the 
entropy production, we can develop a similar argument by considering large, 
negative displacements brought about by negative velocities over short-lived collisions. 



However, in this case it is much simpler to combine the FR of Eq. ( 12 ) and the LDP 



( 36 1 to obtain 

P(As tot /t = -s) ~ e -* s -«V*v^ „ e -* s , s>s*. (38) 

Therefore, similarly to its center, the pdf P(Astot/t = s) satisfies an LDP for 
s < — s* whose dominant exponential scale is t and whose rate function at this scale is 
I(s) — ~s. This result is illustrated in Fig.[ljb) and is consistent with the divergence of 
^Astot W f° r & < ~ 1- Thus, to summarize, we have two LDPs for P(As tot /t — s): one 
at speed t, which correctly characterizes the large deviations of the entropy production 
for s < s* , but not for s > s* since I(s) = in that region; and a second LDP at 
speed y/t, which refines the first one in the region s > s* . 
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6. Conclusions 

We have discussed in this paper the interplay that exists between the fluctuation 
relation (FR) and the large deviation principle (LDP) for a collisional Maxwell-Lorentz 
gas. We have used this model to demonstrate that the FR is more fundamental in 
a sense than the LDP. Focusing on the entropy production of this model, we have 
indeed shown that, although this quantity satisfies an FR, the pdf of this quantity 
involves two different large deviation scalings with time, which implies that an FR 
does not necessarily arise, as often thought, from a uniform LDP decaying everywhere 
exponentially with time. The generality of the FR in this case enabled us to obtain 
corrections to the exponential LDP, thereby showing that the FR can be used in a 
constructive way to obtain large deviations when other approaches, based for example 
on generating functions and the Gartner-Ellis Theorem, fail. 

Although the model studied is greatly simplified compared to real gases, it 



serves as a valuable benchmark to study granular systems 35 36 and conduction 



problems 41 , among other physical phenomena. Moreover, the fact that the non- 
uniform LDP that we find is related to long ballistic accelerations of the gas' particles 
shows that our results should apply to more general collisional models, which might 
include, e.g., a dependence of the mean free time on the relative velocity between 
particles. Recent studies of the hard-sphere case have shown that the distribution 



of times between collisions conserves an exponential tail at large times 42 , so the 
mechanism discussed here is relevant for this case. 



Appendix A. Calculation of the SCGFs 

Appendix A.l. Work 



We can expand from Eq. ( 14 1 the expression of the SCGF of the work W(t) in the 
following way: 



X w (k) = lim \ In [ f>(") (e^ S?=i -)) 

\n=0 ' / 



lim - In 

t— yoo t 



lim - In 

t— s-oo t 

0.(k) - 1 



1 / t 



71=0 



n-2 



E 



1 ftG x (k)\ n 



GAk) 2 ^ nl\ r c 



(A.l) 



where G x (fc) is the generating function of the displacements Xi defined after Eq. |23[ ) . 
In the second line, we have used the statistical independency of the n — 2 displacements 
Xi and the fact that the number n of collisions in a time t is distributed according to 
the Poisson statistics: 



n 



(A.2) 



The displacements are themselves random variables taken from the following pdf: 

r0C dr 



dv t 



-qv 



e-^ T "S[x - (vt + t 2 £/2)}, (A.3) 
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/OO />00 J 

dx e 2kq£x P x {x) = — exp 
-oo JO T c 



q£ 2 T Z k{l + k) 
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(A.4) 



Appendix A. 2. Boundary term 

To find the SCGF of the boundary term, we write it as the sum B — b\ + b n of the two 



terms defined in Eqs. (25), which are asymptotically independent, since they involve 



velocities that are separated by a large number of independent random collisions. As 
a result, we can write 



^kB\ /„kbi 



){e 



kb„ 



(A.5) 



in the limit t — ¥ oo. Of course, the velocity v(0) at the initial time is correlated 
with the velocity v(t^) before the first collision, while the velocity v(t+) after the last 
collision is correlated with the velocity v(t) at the final time. In fact, we have 

«(if) = v(0) + £r 

v(t) = v{t+) + £a, (A.6) 

where r and a are random variables distributed according to the exponential pdf 
P t (t). From this, it is easy to see that the probability distribution of the intervals 
between an arbitrary time, e.g., t = and the following collision t^, is the same as 
that between successive collisions [43]. Using (A.6 1 above, we must then have 



b 1 = ^(v + Sr) 2 + lnP(v) 

b n = -^u 2 - In P(-(u + £a)), 



(A.7) 



where v has the stationary pdf ^ and u has the Gaussian pdf From these pdfs, 
we thus write 



a kb n \ 



dv 



dv 



00 dr 

Tc 



\ 



dr 



e -r/r c _ pk+l^ e 



-t/t c p— 



kq(v+£t) 



£r)e 



-{k+l)qu 2 



where q = m/26 = M/2T. 

At this point, we can use the asymptotic behaviors of -P s tat(iO, given by 

, 3~ bv for v 

"stat(v) 



for v 



-oo 
-oo, 



(A.8) 
(A.9) 

(A.10) 



to find the regions of convergence of the integrals ( A.8 ) and ( A.9 1 . For the first integral, 



we find k € (— oo,0], while for the second one we have k G (—1,0]. These domains 



of convergence determine the domain of convergence of \As tot (k) via Eqs. (A.5) and 
p2l. 
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